Antifreeze proteins and homogeneous nucleation: On the physical determinants impeding ice crystal growth

Antifreeze proteins (AFPs) are biopolymers capable of interfering with ice growth. Their antifreeze action is commonly understood considering that the AFPs, by pinning the ice surface, force the crystal–liquid interface to bend forming an ice meniscus, causing an increase in the surface free energy and resulting in a decrease in the freezing point ΔT max. Here, we present an extensive computational study for a model protein adsorbed on a TIP4P/Ice crystal, computing ΔT max as a function of the average distance d between AFPs, with simulations spanning over 1 μs. First, we show that the lower the d, the larger the ΔT max. Then, we find that the water–ice–protein contact angle along the line ΔT max(d) is always larger than 0○ , and we provide a theoretical interpretation. We compute the curvature radius of the stable solid–liquid interface at a given supercooling ΔT ≤ ΔT max, connecting it with the critical ice nucleus at ΔT. Finally, we discuss the antifreeze capability of AFPs in terms of the protein–water and protein–ice interactions. Our findings establish a unified description of the AFPs in the contest of homogeneous ice nucleation, elucidating key aspects of the antifreeze mechanisms and paving the way for the design of novel ice-controlling materials

Estudio Integral sobre la Evaluación del Recurso Eólico en Entornos Urbanos. Estación Anemométrica Adaptable

Como consecuencia de un intenso proceso de trabajo, coordinado con Investigadores del CIEMAT (Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas), se pretende comunicar el estado del trabajo de Investigación relacionado con la Evaluación del Potencial Eólico en entornos urbanos, con vistas a la posible generación de energía eléctrica limpia

Solubility of Methane in Water: Some Useful Results for Hydrate Nucleation

In this paper, the solubility of methane in water along the 400 bar isobar is determined by computer simulations using the TIP4P/Ice force field for water and a simple LJ model for methane. In particular, the solubility of methane in water when in contact with the gas phase and the solubility of methane in water when in contact with the hydrate has been determined. The solubility of methane in a gas–liquid system decreases as temperature increases. The solubility of methane in a hydrate–liquid system increases with temperature. The two curves intersect at a certain temperature that determines the triple point T3 at a certain pressure. We also determined T3 by the three-phase direct coexistence method. The results of both methods agree, and we suggest 295(2) K as the value of T3 for this system. We also analyzed the impact of curvature on the solubility of methane in water. We found that the presence of curvature increases the solubility in both the gas–liquid and hydrate–liquid systems. The change in chemical potential for the formation of hydrate is evaluated along the isobar using two different thermodynamic routes, obtaining good agreement between them. It is shown that the driving force for hydrate nucleation under experimental conditions is higher than that for the formation of pure ice when compared at the same supercooling. We also show that supersaturation (i.e., concentrations above those of the planar interface) increases the driving force for nucleation dramatically. The effect of bubbles can be equivalent to that of an additional supercooling of about 20 K. Having highly supersaturated homogeneous solutions makes possible the spontaneous formation of the hydrate at temperatures as high as 285 K (i.e., 10K below T3). The crucial role of the concentration of methane for hydrate formation is clearly revealed. Nucleation of the hydrate can be either impossible or easy and fast depending on the concentration of methane which seems to play the leading role in the understanding of the kinetics of hydrate formation

Melting points of water models: Current situation

By using the direct coexistence method, we have calculated the melting points of ice Ih at normal pressure for three recently proposed water models, namely, TIP3P-FB, TIP4P-FB, and TIP4P-D. We obtained Tm = 216 K for TIP3P-FB, Tm = 242 K for TIP4P-FB, and Tm = 247 K for TIP4P-D. We revisited the melting point of TIP4P/2005 and TIP5P obtaining Tm = 250 and 274 K, respectively. We summarize the current situation of the melting point of ice Ih for a number of water models and conclude that no model is yet able to simultaneously reproduce the melting temperature of ice Ih and the temperature of the maximum in density at room pressure. This probably points toward our both still incomplete knowledge of the potential energy surface of water and the necessity of incorporating nuclear quantum effects to describe both properties simultaneously.Depto. de Química FísicaFac. de Ciencias QuímicasTRUEpu

Interfacial Free Energy and Tolman Length of Curved Liquid-Solid Interfaces from Equilibrium Studies

In this work, we study by means of simulations of hard spheres the equilibrium between a spherical solid cluster and the fluid. In the NVT ensemble we observe stable/metastable clusters of the solid phase in equilibrium with the fluid, representing configurations that are global/local minima of the Helmholtz free energy. Then, we run NpT simulations of the equilibrated system at the average pressure of the NVT run and observe that the clusters are critical because they grow/shrink with a probability of 1/2. Therefore, a crystal cluster equilibrated in the NVT ensemble corresponds to a Gibbs free energy máximum where the nucleus is in unstable equilibrium with the surrounding fluid, in accordance with what has been recently shown for vapor bubbles in equilibrium with the liquid. Then, within the Seeding framework, we use Classical Nucleation Theory to obtain both the interfacial free energy γ and the nucleation rate. The latter is in very good agreement with independent estimates using techniques that do not rely on Classical Nucleation Theory when the mislabeling criterion is used to identify the molecules of the solid cluster. We therefore argue that the radius obtained from the mislabeling criterion provides a good approximation for the radius of tension, R_s . We obtain an estimate of the Tolman length by extrapolating the difference between R e (the Gibbs dividing surface) and R s to infinite radius. We show that such definition of the Tolman length coincides with that obtained by fitting γ versus 1/R_s to a straight line as recently applied to hard spheres

“In Silico” Seawater

Many important processes affecting the earth’s climate are determined by the physical properties of seawater. In addition, desalination of seawater is a significant source of drinking water for the human population living in coastal areas. Since the physical properties of seawater governing these processes depend on the molecular interactions among its components, a deeper knowledge of seawater at the molecular level would contribute to a better understanding of these phenomena. However, in strong contrast with the situation in other areas such as biomolecules or materials science, molecular simulation studies reporting the physical properties of seawater are currently lacking. This is probably due to the usual perception of the seawater composition being too complex to approach. This point of view ignores the fact that physical properties of seawater are dependent on a single parameter representing the composition, namely the salinity. This is because the relative proportions of any two major constituents of seasalt are always the same. Another obstacle to performing molecular simulations of seawater could have been the unavailability of a satisfactory force field representing the interactions between water molecules and dissolved substances. However, this drawback has recently been overcome with the proposal of the Madrid-2019 force field. In this work we show for the first time that molecular simulation of seawater is feasible. We have performed molecular dynamics simulations of a system, the composition of which is close to the average composition of standard seawater and with the molecular interactions given by the Madrid-2019 force field. In this way we are able to provide quantitative or semiquantitative predictions for a number of relevant physical properties of seawater for temperatures and salinities from the oceanographic range to those relevant to desalination processes. The computed magnitudes include static (density), dynamical (viscosity and diffusion coefficients), structural (ionic hydration, ion−ion distribution functions), and interfacial (surface tension) properties

A simulation study of homogeneous ice nucleation in supercooled salty water

We use computer simulations to investigate the effect of salt on homogeneous ice nucleation. The melting point of the employed solution model was obtained both by direct coexistence simulations and by thermodynamic integration from previous calculations of the water chemical potential. Using a seeding approach, in which we simulate ice seeds embedded in a supercooled aqueous solution, we compute the nucleation rate as a function of temperature for a 1.85 NaCl mol per water kilogram solution at 1 bar. To improve the accuracy and reliability of our calculations, we combine seeding with the direct computation of the ice-solution interfacial free energy at coexistence using the Mold Integration method. We compare the results with previous simulation work on pure water to understand the effect caused by the solute. The model captures the experimental trend that the nucleation rate at a given supercooling decreases when adding salt. Despite the fact that the thermodynamic driving force for ice nucleation is higher for salty water for a given supercooling, the nucleation rate slows down with salt due to a significant increase of the ice-fluid interfacial free energy. The salty water model predicts an ice nucleation rate that is in good agreement with experimental measurements, bringing confidence in the predictive ability of the model. We expect that the combination of state-of-the-art simulation methods here employed to study ice nucleation from solution will be of much use in forthcoming numerical investigations of crystallization in mixtures. Published by AIP Publishing

Building a Hofmeister-like series for the maximum in density temperature of aqueous electrolyte solutions.

The temperature of the maximum in density (TMD) at room pressure is experimentally evaluated for aqueous solutions of a set of halides containing F−, Cl−, Br−, I−, Li+, Na+, K+, Rb+, Cs+ and Mg2+ at a 1 m concentration. The measurements were performed by monitoring the density- temperature profiles and tracking the temperature-dependent position of the meniscus, in a capillary glass tube. Adding salts diminishes the TMD of the solutions with respect to pure water, being the magnitude of the change dependent on the nature of the electrolyte. The experimental values of the shift in the TMD can be split into individual ion contributions. From this information we were able to establish a rank of ions (i.e. a Hofmeister-like series) according to their efficiency in shifting down the TMD. The experimental results are also compared to simulation values obtained via Molecular Dynamics using the Madrid-2019 force field that assigns non-integer charges for the ions and is parametrized for the TIP4P/2005 water model. Finally, since the TMD is a fingerprint property of water, we will discuss the impact of ions on this maximum in relation with the way different ions modify the structure of water